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Abstract 


In  many  chemical  engineering  control  problems  the 
amount  of  movement  of  the  manipulated  variable  is  as  impor¬ 
tant  as  that  of  the  controlled  variable.  In  this  study  such 
a  case  was  examined.  The  system  investigated  was  one  in  which 
the  controlled  variable  had  to  only  stay  within  limits.  The 
distance  between  the  limits  was  then  used  to  lower  the  fluctua¬ 
tions,  root-mean-square,  of  the  manipulated  variable.  A  study 
of  nonlinear  controllers  was  carried  out  to  see  if  the  root- 
mean-square  could  be  lowered  while  still  preserving  a  rapid 
and  highly  damped  start-up  response  or  mathematically  a  low 
value  of  the  integral  with  respect  to  time  of  the  absolute 
value  of  the  manipulated  variable  times  time.  Using  a 
statistical  describing  function  approach  for  analytical  root- 
mean-square  calculations  and  an  analogue  computer,  different 
continuous  nonlinearities  were  investigated.  The  nonlinearity 
with  the  best  root-mean-square  and  above  integral  values,  in 
the  author's  opinion,  is  one  in  which  a  simply  proportional 
controller  has  its  input  changed  from  a  linear  error  to 
j(error)  I  x  (error)  and  the  gain  of  the  proportional  controller 


itself  is  lowered. 
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Nomenclature 

2 

area  (length) 
error, (length) 

mean-square  of  error  (length)2 
root-mean-square  of  error  (length) 
liquid  level  of  tank  (length) 
mean-square  of  liquid  level  (length)2 
root-mean-square  of  liquid  level  (length) 
statistical  describing  function  (length/length) 
gain 

1/area  (1/length2) 
probability  density  (1/length) 

3 

flow  out  (manipulated  variable)  (length)  /sec. 
flow  in  (load  variable)  (length)  /sec. 
reference  (length) 
time  (seconds) 

any  fluctuating  variable  (length) 

2 

power  spectral  density  of  v  (length)  (seconds) 

2 

power  spectral  density  of  x  (length)  (seconds) 


Symbols  for  Laplace  Transforms  of  Variables 
E  Laplace  transform  of  e 

h(s)  Laplace  transform  of  h 

Q  Laplace  transform  of  qQ 

Qs  Laplace  transform  of  q^ 
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Introduction 

The  theory  of  linear  control  is  well  developed.  If  a 
system  and  its  controller  are  linear,  one  needs  only  adjust 
the  proportional,  reset,  or  rate  constants  to  obtain  whatever 
criterion  was  specified  such  as  minimum  area,  stability,  etc. 
Once  a  criterion  has  been  found  it  can  be  transferred  to 
related  systems  with  ease.  However,  few  systems  in  practice 
are  truly  linear.  Also,  in  many  cases,  it  is  possible  to  de¬ 
sign  a  nonlinear  controller  that  will  yield  a  better  response 
than  a  linear  controller.  As  a  result,  nonlinear  control 
theory  has  been  in  much  demand. 

Using  nonlinear  control  theory,  it  is  generally 
possible  to  design  a  better  controller  for  two  reasons.  First, 
instead  of  replacing  the  nonlinearity  with  an  approximate 
linear  expression  or  assuming  it  does  not  exist,  the  non¬ 
linearity  can  be  studied  for  its  true  effect  on  the  closed 
loop.  Second,  one  can  then  design  the  nonlinear  controller 
to  compensate  or  use  nonlinearities  to  obtain  the  best  possible 
response.  Thus,  it  might  be  said  that  nonlinear  control  theory 
gives  one  the  power  to  tailor  each  controller  more  closely  to 
a  system's  particular  needs. 

Nonlinear  control  theory  has  been  slow  in  general 
acceptance  because  of  the  difficulty  normally  encountered  in 
obtaining  an  analytical  solution  to  a  control  problem.  For  up 
to  the  present  time,  there  is  no  general  analytical  approach 
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to  all  nonlinear  problems.  Instead,  there  are  certain  types 
of  approximations,  such  as  the  describing  function  approach  to 
stability  analysis  or  statistical  describing  function  for  rms 
(root— mean-square)  calculations  as  used  in  the  present  study; 
and  then  there  are  special  approaches,  such  as  the  phase  plane 
technique  which  is  used  to  check  the  response  of  second  order 
systems.  If  the  feedback  loops  and  nonlinearities  become 
numerous,  an  analytical  solution  quickly  becomes  impractical 
if  not  impossible.  When  this  occurs  a  solution  is  obtained  by 
using  the  computer,  either  digital  or  analogue. 

With  the  above  thoughts  in  mind,  this  study  designs  a 
nonlinear  controller.  The  problem  is  such  that  it  is  desired 
to  keep  the  controlled  variable  within  limits,  only,  while 
keeping  the  manipulated  variable  from  fluctuating  as  much  as 
possible.  Also,  the  controlled  variable  must  not  feel  large 
upsets  for  long.  To  explain  this  more  fully,  consider  the 
reboiler  on  a  distillation  column.  Here,  it  is  desired  to  keep 
the  liquid  level  in  the  reboiler  above  the  steam  tubes  and  a 
given  distance  below  the  bottom  plate.  At  the  same  time 
variations  in  the  outflow  must  be  minimized  so  as  not  to  dis¬ 
turb  downstream  processes.  Also,  sudden  load  disturbances  such 
as  occur  during  start-up  should  damp  out  quickly.  This  problem 
will  be  investigated  throughout  the  remainder  of  this  study. 

One  must  now  translate  the  above  statements  of  desired 
response  into  definite  mathematical  criteria.  After  the 
criteria  are  chosen  an  order  of  priority  must  then  be  established. 
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Both  of  these  decisions  are  dependent  on  the  designer's 
judgement.  The  criteria  that  guided  the  outcome  of  this 
study  are  listed  below  in  order  of  priority  as  follows: 

1)  Stability 

2)  The  lowest  value  of  ITAMff  rqQ  tdt)  for  step 
changes  in  the  liquid  level  (qQ  =  manipulated  variable, 
t  *  time) 

3)  The  lowest  rms  (root-mean-square)  value  for  the 
manipulated  variable  for  a  random  Gaussian  load  change. 

Continuous  functions  are  studied  to  meet  the  above 
demands . 
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Literature  Review 

The  following  short  summaries  and  reviews  represent 
the  literature  covering  the  subject  of  this  study.  The  first 
paper  is  one  showing  how  this  problem  can  be  handled  in  the 
design  stage.  The  second  paper  describes  a  commercial  con¬ 
troller  that  is  used  to  smooth  out  flows  from  tanks.  In  the 
last  part  of  this  summary  and  review  a  discussion  is  presented 
of  other  controllers  and  theory  in  general. 

R.  Jackson  (5)  suggests  that  the  object  of  all  con¬ 
trollers  is  to  maintain  some  measured  quantity  y(t)  close  to 
Y,  in  spite  of  a  disturbance  i(t).  He  then  mentions  that  the 
following  work  arose  because  it  was  felt  that  a  lack  of  infor¬ 
mation  at  the  design  stage  should  not  lead  to  complete  im¬ 
potence.  Rather,  the  system  should  be  over  designed  based  on 
a  safety  factor  that  is  related  to  the  designer's  ignorance 
of  the  system. 

Jackson  then  points  out  that  at  the  design  stage  the 
maximum  values  of  quantities  and  rates  are  generally  known; 
and  he  then  goes  on  to  state  the  basic  mathematical  problem  as 
follows:  "In  a  stable  linear  filter  the  output  y(t)  at  time  t 
is  related  to  values  of  the  input  i(t)  at  all  previous  times 
by  an  equation  of  form 

00 

y(t)  =  {  W (u)  i(t  -  u)  du  (1) 

0 

where  the  weighting  function  W(u)  characterizes  the  particular 
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system  considered  and  the  stability  conditions  demand  that 

CO 

/  |  W(u) I du  should  exist,  It  is  assumed  that  i(t)  is  uniformly 
bounded  and  that  the  origin  for  i(t)  is  chosen  so  that 

U(t)  |  *  im  (all  t)  (2) 

It  is  also  assumed  that  the  rate  of  change  of  i (t)  is  uniformly 
bounded  or,  slightly  more  generally,  that 

I  i  (t1)  •-  i  (t„)  1 

- £ .  ^  (all  tx,t2)  (3) 

~  L2  ; 

The  problem  is  then  to  find  the  corresponding  smallest  bound 
ym  such  that  !  Y  ( t )  |  4'  for  all  t."  (  )  A  procedure  is  then 
developed  for  constructing  that  i (t)  which  makes  y ,  at  some 
specified  time,  as  large  as  possible.  The  value  of  y(=  y  ) 
can  then  be  obtained  from  equation  (1) .  The  solution  to  this 
problem  can  be  found  when  the  performance  is  specified  by  limits 
±y  and  when  the  information  available  about  the  disturbance 
is  of  the  form  specified  by  equations  (2)  and  (3) . 

Jackson  then  proceeds  to  develop  the  mathematical 
theory  for  solution  in  a  rigorous  manner.  Having  formally 
proved  his  method,  he  then  applies  his  procedure  to  the 
practical  application  of  level  control  systems.  He  divides 
the  types  of  level  control  into  three  classes  called  flow 
smoothing,  level  control,  and  given  valve  motor. 
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Figure  1 

The  flow  smoothing  system,  as  described  by  Jackson, 
makes  use  of  the  available  free  space  in  the  vessel  to  smooth 
out  fluctuations  ir  the  out  flow  f ^  (see  Figure  1) .  The  aim 
here  is  to  vary  as  slowly  as  possible  while  still  checking 
the  level  variations  before  they  overflow  or  empty  the  tank*, 
Jackson  states  the  object  of  design  is  to  predict: 

1,  The  minimum  size  of  vessel  which  can  be  used  if 
the  desired  smoothness  of  the  controlled  flow  is  to  be  obtained 
without  any  possibility  of  the  vessel  overflowing  or  emptying . 

2 „  Size  and  speed  requirements  for  the  control  valves 
regulating  the  flow 

With  the  level  control  system,  it.  is  desired  to  keep 
the  liquid  level  between  narrow  limits •  Here  the  vessel  size 
is  assumed  given  and  the  designer  must  specify  a  valve  of 
adequate  speed  and  size.  The  third  control  class,  given  valve 
motor,  is  the  reverse  problem  of  the  level  control  system. 

All  the  above  control  problems  are  basically  the  same, 
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Thus,  assuming  proportional  control,  Jackson  develops  graphs 
so  that  the  solution  to  all  the  above  problems  can  be  obtained 
without  the  considerable  amount  of  work  that  would  otherwise 
be  needed. 

This  short  summary  of  Jackson's  work  thus  shows  how 
the  problem  of  this  study  could  be  handled  in  the  design 
stage . 

Existing  Controller 

Foxboro  Company  has  developed  a  nonlinear  tool  which 
makes  a  dead  zone  around  the  reference  to  help  smooth  out  the 
flow(  ) .  This  device  is  used  in  conjunction  with  a  proportional 
or  proportional-reset  controller. 

There  is  no  question  that  a  controller  with  a  dead 
zone  will  most  efficiently  lower  the  rms  of  the  manipulated 
variable  if  the  dead  zone  is  made  relatively  large.  Should 
there  be  an  upset,  however,  would  this  controller  be  able  to 
handle  it?  Can  it  handle  start-up?  The  answer  to  these 
questions  and  the  related  ones  of  relative  stability  and  damping 
properties  can  best  be  answered  by  remembering  the  phase  plane 
diagram  of  a  second  order  system  with  a  dead  zone.  On  the  phase 
plane  the  dead  zone  is  seen  to  be  a  hinderance  to  damping  for 
step  inputs.  The  larger  the  dead  zone  the  greater  the 
hinderance.  Since  it  is  desired  to  have  both  the  ability  to 
damp  out  sudden  changes  and  to  lower  the  rms,  this  type  of 
controller  can  be  improved  upon. 
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The  great  advantage  of  this  controller  is  its  versatility 
for  it  is  used  in  many  applications  besides  the  above-mentioned 
one. 

General  Theory 

O.J.Me  Smith  has  considered  many  control  devices.  One 
which  was  of  interest  was  posicast  control (7).  This  is  because 
posicast  control  makes  the  manipulated  variable  react  swiftly 
to  step  changes  in  the  controlled  variable  and  does  not  let  the 
system  fluctuate  in  the  ideal  use.  However,  this  control  is 
of  little  value  here  for  two  reasons.  One  is  that  it  is  limited 
to  underdamped  systems.  Two,  the  rms  of  the  manipulated 
variable  would  be  much  too  high  if  it  responded  to  each 
fluctuation  of  the  controlled  variable. 

O.J.M.  Smith  has  covered  much  nonlinear  control  theory 
in  general (7).  In  one  section,  he  shows  how  to  use  the  phase 
plane  for  the  design  of  servo  control.  While  illustrating  his 
method  with  an  example  problem,  he  uses  the  |  (variable  lx  (variable) 
expression  used  in  this  study  as  eje|. 

Brown  and  Nilsson (  )  give  an  excellent  introduction 
into  the  use  of  random  signals  as  applied  to  linear  systems. 

Thaler  and  Pastel (  )  advance  the  discussion  of  the  use  of 
random  signals  to  nonlinear  systems;  and  it  is  their  methods 
together  with  those  presented  by  Gibson (  )  which  are  used  in 
this  study. 
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Theory 

In  the  theoretical  development  a  mathematical  model 
for  the  process  dynamics  is  first  obtained.  A  control  loop  is 
then  established  and  the  reasons  for  a  given  filter  and  for 
discarding  certain  nonlinearities  are  given.  A  stability 
analysis  is  then  carried  out  on  both  the  linear  and  nonlinear 
systems 0  Finally,  a  short  mathematical  development  is  given 
of  Thaler  and  Pastel's  (  )  method  of  obtaining  rms  values  in 
linear  and  nonlinear  systems  that  have  been  excited  by  random 
load  fluctuations.  The  author  then  extends  the  method  to 
study  the  response  to  random  disturbances  of  the  manipulated 
variable  as  well  as  the  controlled  variable. 

The  mathematical  model  is  developed  using  Figure  2. 


Figure  2 
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Writing  a  material  balance 


dh 

A  —  =  q 

dt  1 


The 


A 

t 

h 


Laplace 


~  tank  area 
=  time 

=  liquid  level 
=  flow  in 
=  flow  out 

transform  of  this  equation  may  be  written 


A  S  h  (S )  =  L  [qj 


L  M 

h  (S)  =  - — 

A  S 


&o] 


provided  that  h  |o+j  is  assumed  equal  to  zero.  A  capacitance 
element  thus  describes  the  dynamic  response  of  the  process. 


The  control  loop  is  now  closed  giving  the  block  diagram  of 


Figure  3* 


Figure  3 
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The  selection  of  a  filter  is  arbitrary*  However,  a  first 
order  filter  with  a  time  constant  of  0,1  of  A  will  remove  un¬ 
wanted  high  frequencies  in  the  feedback  signal;  yet,  it  will 
not.  seriously  interfere  with  the  response  of  the  system.  Thus,  the 

transfer  function  of  the  filter  is  =— where  T_  =  0.1A, 

1  Hr  1  r 

F 

It  is  intended  to  operate  the  system  within  its  flow 
limits;  therefore,  the  nonlinearity,  saturation  of  the  valve 
output,  is  dropped  from  the  model.  The  semi-log  valve  character¬ 
istic  is  also  neglected.  This  is  because  the  nonlinearity  is 
minimized  by  the  pressure  drop  relationship  between  the  valve 
and  the  pipe  as  the  flow  changes.  Other  nonlinearities,  such 
as,  stiction  in  the  valves  are  neglected  because  of  the  very 
small  effect  they  have  on  the  process.  For  the  reasons  given 
above  the  block  diagram  of  Figure  3  can  simplified  to  that  of 
Figure  4, 


Figure  4 
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The  block  diagram  of  Figure  4  represents  the  system  studied 
in  this  work.  At  this  stage,  the  author  would  like  to  emphasize 
the  fact  that  the  controlled  variable  h  has  only  to  stay  within 
limits.  It  is  then  intended  to  use  the  region  between  these 
limits  to  damp  out  fluctuations  in  the  manipulated  variable  qQ. 
Also,  it  is  desired  that  the  ITAM  ({  |qQ|tdt)  be  small  for  large 
steps  because  of  start-up  and  possible  large  upsets. 

Stability  Studies 

If  linear  control  is  used  on  the  system  of  Figure  4,  it 
is  found,  by  computer,  that  a  simple  proportional  element  with 
a  gain  of  four  gives  the  lowest  ITAM  value.  For  these  con¬ 
ditions,  the  block  diagram  is  that  of  Figure  5. 


Figure  5 

The  Root  Locus  plot  of  the  characteristic  equation  of 


Figure  5  is  shown  in  Figure  6. 
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Figure  6 

Because  no  roots  lie  in  the  right-half -plane  the  linear  system 
is  stable  for  all  gains.  Note,  however,  that  as  the  gain 
increases  the  system  becomes  more  oscillatory. 

A  study  of  the  stability  of  the  nonlinear  functions 
e  |  e |  and  e^  that  proved  to  be  of  interest  is  now  carried  out 
using  the  describing  function  approach(9).  This  method  is  an 
approximation  which  is  based  on  the  assumptions  that: 

1)  The  input  to  the  nonlinear  component  is  a  pure 
sine  wave. 

2)  The  output  curve  contains  no  zero  frequency  com¬ 
ponent  and  no  subharmonic  terms,  and  all  higher  harmonic 
terms  may  be  discarded. 

A  procedure  suggested  when  using  this  method  is  as 
follows : 

1)  Isolate  the  nonlinear  component  in  the  block  diagram. 

2)  Develop  a  Fourier  series  representation  of  the  non¬ 
linearity  based  on  the  above  assumptions. 
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3)  The  describing  function  is  then  the  ratio  of  the 
magnitude  of  the  fundamental  term  in  the  Fourier  series  for 
the  output  wave  to  the  magnitude  of  the  input  sinusoid  at  a 
phase  angle  which  is  the  angle  between  the  two  sine  waves. 

4)  The  product  of  the  describing  function  and  the 
linear  characteristic  equation  is  then  studied  in  the  normal 
fashion  for  stability.  (When  the  describing  function  is 
frequency  dependent,  a  new  Root  Locus  plot  etc.  must  be  con¬ 
structed  for  each  frequency.) 

In  mathematical  notation,  if  the  input  wave  is 
designated  by 


Input  =  E  sin  wt 

and  the  fundamental  frequency  term  in  the  Fourier  series  for 
the  output  is  designated  by 


Fundamental  output  term  =  F(E,w)  sin  [wt  +  <J>(E,wf]  then 
the  describing  function  for  the  nonlinear  component  is  defined 
to  be 


Describing  function  =  GD(E,jw) 


F(E,w) 


<P  (E  ,w) 


The  above  method  is  now  applied  to  the  system  of  Figure 
4  which  has,  first,  a  nonlinear  controller  with  a  K  e|e| 

3 

characteristic  and,  second,  a  controller  with  a  K  e  character- 

2 

istic.  Letting  g(e)  =  e|e|  or  e  and  g(  )  equal  the  mathe- 
matical  operations  | variable  |x  variable  and  (variable)  ,  the 
block  diagram  becomes  that  of  Figure  7 . 
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H 

1 

h 

TpS  +  1 

Figure  7 


The  nonlinearities  are  already  isolated  in  this  case e  There¬ 
fore,  the  Fourier  series  and  describing  functions  are  easily 
obtained  for  the  two  nonlinearities  as  outlined  in  the  above 
procedure  a 

a)  Calculations  for  the  describing  function  of  the 
nonlinearity  e|e| 


e  =  E  sin  wt 
2 

e | e |  =  E  sin  wtjsin  wt! 

The  two  constants  for  the  fundamental  frequency  in  the  Fourier 
series  are: 


2  it 


A, 


t  2  2 

J  E  sin  wt  cos  wt  dwt 


TT  0 

0 


Note:  sin  wt  -  sin  wt  from  0  to 


B. 


2  Tf  ^  o 

—  J  Ez  sin"  wt  dwt 


it  o 
4E2 


0*425  E' 


3  IT 


s  :  9  l&ixeaaUnon 
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Therefore 

0.425E2 

G  (E,jw)  =  - -  =  0.425  E 

E 


b)  Calculations  for  the  describing  function  of  the 
3 


e  =  E  sin  wt 


3  3  3 

e  =  E“  sin  wt 


Again  the  Fourier  series  constants  are: 


IT 


3  3 

—  /  E  sin^  wt  cos  wt  dwt 

Tf  0 


0 


B. 


2  TT 

t  3  «  4 

—  j  E  sm  wt  dwt 

71  0 


Therefore 


-  E 
4 


E' 


GD(E,jw)  -  0.75 


=  0.75  E 


The  calculations  done  in  a)  and  b)  above  show  that  the 

i  i  3 

describing  function  for  both  e|e|  and  e  are  only  amplitude 
dependent o  Thus,  on  applying  step  4  of  the  suggested  pro¬ 
cedure,  it  can  be  seen  that  both  nonlinear  systems  are  stable 
because  the  linear  system  is  stable  for  all  gains. 
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Root -Me an -Square  (rms)  Calculations 

The  theory  is  first  presented  for  the  mean-square  (ms) 
calculations  in  the  frequency  domain  for  any  random  variable 
v.  Second,  a  method  is  discussed  that  uses  a  statistical 
describing  function  to  deal  with  nonlinearities  when  the  system 
is  under  random  load  fluctuations .  Third,  the  theory  behind 
the  statistical  describing  function  is  statedc  Fourth,  the 


above  methods  of  calculating  the  rms  are  applied  to  the  con¬ 
trolled  variable  h.  Fifth,  the  author  extends  the  rms  calcula¬ 
tions  to  the  manipulated  variable  q  .  Sixth,  the  author  pre¬ 
dicts  approximately  the  effect  given  nonlinearities  will  have  on 
both  h  and  q  rms  values. 


A  random  variable  which  has  been  described  in  terms  of 
its  time  variation  may  also  be  described  in  terms  of  its 
frequency  spectrum.  The  following  thus  develops  a  relation¬ 
ship  between  the  mean-square  (ms)  value  of  a  random  variable 
and  its  frequency  content. 


II 


A  convenient  starting  point  is  the  time  equation  for 


the  ms  value  which  is 


v 


ms 


lim 

T-*°o  2T  -T 


(1) 


v 


variable,  T  -  time  limit,  t  =  time 


There  are  obvious  mathematical  advantages  available  if  one  or 
both  of  the  variable  functions  v(t)  in  equation  (1)  can  be 
replaced  by  their  Fourier  transforms.  Unfortunately,  this 
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cannot  be  done  directly  since  the  Fourier  integral  of  a  random 
time  function  does  not  converge .  However,  if  a  new  function 
vT(t)  is  defined,  the  desired  manipulation  may  be  made.  The 
function  vT(t)  is  defined  to  exist  only  in  the  time  interval 
-T  <  t  <  T  ,  being  exactly  the  same  as  v(t)  during  this  inter¬ 
val  but  being  zero  outside  the  interval.  With  this  definition 
(1)  may  be  rewritten 


1  T 


v 


ms 


lim  —  J 

T*oo  2T  -T 


(2) 


The  Fourier  transform  of  the  variable  v^Ct)  exists 


T 


1 


Substituting  (3)  in  (2)  for  only  one  factor  of  vT(t) 


Changing  the  order  of  integration 


v 


ms 


IT  .  . 

V  (w)  C  lim  —  /  vT(t)  e 3  dt]  dw  (5) 


•p-^oo  2  T  — T 


but 


CO 


V  (w)  =  J  vT(t)  e  jwt  dt 

1  no  A 


(6) 


Therefore,  in  equation  (  )  it  is  seen  that 


T 

lim  / 


dt 


(7) 


. 
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Thus  equation  (5)  may  be  rewritten 


v 


ms 


1  00  1 

—  J  lim  —  [Vrn(w)  V  (-W)H  dw  (8) 
2v  -«*  T-*<*  2T  1  1 


From  equation  (8)  the  power  spectral  density  is  defined  as 

1 


*v  (w) 


lira 

•Phm» 


2T 


L  VT(w)  Vwt-*Q 


(9) 


It  is  readily  seen  that  equation  (8)  approaches  equation  (2) 
as  T-*00 .  Hence,  the  desired  relationship  between  the  mean- 
square  (ms)  value  of  the  random  variable  and  its  frequency 
content  is  (4)" 


v 


ms 


1  <* 

- —  /  (w)  dw 

2  w  -<* 


(10) 


The  power  spectral  density  function  equation  (9)  is 
simply  the  Fourier  transform  of  the  autocorrelation  function (5)* 
The  autocorrelation  function  is  the  equation  describing  the 
statistical  behavior  of  the  variable  v . 

If  x(t)  is  a  time-stationary  random  function  with  a 

f 

known  autocorrelation  function  (or  power  spectral  density)  and 
the  poles  of  Y(S),  in  Figure  8,  lie  in  the  left-half  plane 
then (5) 

2 


*v(w)  =  *x<w>  1 Y ( jw) | 


where  <j>  (w)  -  power  spectral  density  of  variable  x 


(11) 


x  (t) 

Y  (s') 

v  ( t ) 

Figure  8 
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(The  above  statement  is  proven  in  reference  i#)  This  completes 
the  development  required  to  obtain  ms  and  thus  rms  values  for 
linear  systems. 

The  ms  values  for  nonlinear  systems  are  calculated 

in  a  similar  manner.  The  nonlinearity  is  simply  replaced  by 

the  term  called  the  statistical  describing  function  and 

the  calculations  are  then  carried  out  as  explained  above.  The 

result  for  v  then  contains  the  t  o  m  \\  .  h  is  then 

ms  o  o 

eliminated  by  obtaining  it  as  a  function  of  v  and  solving  the 
equation  for  v  either  analytically  or  graphically. 

It  now  remains  to  be  shown  that  the  statistical  de¬ 
scribing  function  can  be  calculated  as  a  function  of  v. 
First,  it  is  assumed  that  the  load  fluctuations  and  all  other 
variables  in  the  control  loop  including  v,  have  a  Gaussian 
probability  density.  (The  validity  of  this  assumption  is  dis¬ 
cussed  in  reference  9.)  In  this  study  the  average  value  of  all 
variables  in  the  control  loop  are  taken  as  zero.  Therefore, 
the  probability  function  has  the  shape  shown  in  Figure  9 


Figure  9 
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and  is  expressed  mathematically  as 


Probability  density  of  v  =  P(v)  = 

V2 

2  a2 
e 

-  (12) 

a  /2  it 

where 


o  =  standard  deviation  of  v 

When  the  probability  density  is  known, 

the  average  value 

of  the  random  signal  is  by  definition 


00 

vav  =  J  v  P  (v)  dv 

—  00 

(14) 

and  similarly,  the  mean-square  (ms)  value  is 


00 

2 

Vffl(,  =  J  v  P(v)  dv 

ms  ' 

—  00 

(15) 

2 

Also,  if  instead  of  v  in  Equation  (15) ,  we  have  some  other 


function  f(v);  then  its  average  value  f(v)  is 

given  by: 

00 

f(v)  =  J  f(v)  P(v)  dv 

—  00 

(16) 
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Equation  (16)  now  makes  it  possible  to  calculate  the  statistical 
describing  function  Ho„  If  for  the  nonlinear  element,  the 
function  f(v)  is  set  equal  to  the  square  of  the  output  and  vms 
is  the  square  of  the  input  Ho  becomes: 


H 


o 


f  (V) 


(17) 


\  v 
V  ms 

Thus,  the  statistical  describing  function  is  equal  to  the 
square  root  of  the  average  mean  square  of  the  output  over 

r 

input  for  any  nonlinear  element. 

The  above  theory  is  now  applied  to  the  systems  shown 
in  Figures  5  and  7.  When  talking  of  a  statistical  describing 
function,  the  linear  system  in  Figure  5  merely  becomes  a  special 
case  of  that  shown  in  Figure  7.  Therefore,  it  is  necessary 
to  consider  only  the  nonlinear  system  of  Figure  7„ 

Let 


H 

o 

E(S) 

hf  ( s) 

h  (s) 
QCs) 
QD(s) 


=  statistical  describing  function  of  nonlinearities 
«  L  [e  ] 

=  L  [H  ] 

=  L  i>  ] 

=  L  Lqol 

=  L  0?i  1 
1 

A 


The  new  block  diagram  for  the  nonlinear  system  is  shown  in 

Figure  10 „ 
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Figure  10 

Qd(S)  is  the  Laplace  transform  of  the  time  equation 

(intercorrelation)  describing  the  Gaussian  random  movement  of 

the  load  variable.  For  the  mathematical  model  on  the  computer, 

Q^(S)  was  obtained  from  a  Gaussian  random  acnerator  which 

passed  its  output  through  a  low  frequency  filter  with  the  La- 

S 

place  transform  ( - - - )  .  The  Gaussian  random  generator  has 

10S  +  1 

a  high  requency  filter  with  the  Laplace  transform  ( - — - - ) 

0.0159S  +  1 

However,  this  filter  can  be  neglected  for  simplicity.  This  is 

because  the  filter  in  the  feedback  of  the  control  loop,  for 

practical  purposes,  eliminates  all  high  frequencies.  Its 

Laplace  transform  is  ( - - - )  on  the  time  scaled  computer 

0.5S  +  1 

model.  Thus,  before  the  filter  on  the  random  generator  has 

any  effect,  the  high  frequencies  have  been  reduced  at  least 

one  order  of  magnitude  by  the  feedback  loop  filter. 

The  Gaussian  random  generator  output  can  be  represented 

nw 

by  the  Laplace  transform!  — -  where: 

S  +  w 

o 


;  vr-)  .0 

*  v<  I  i>  r:  r<  ,  :U;  a.  brBBBq 


! 


. 


:  £■ 


£•:' .  ■;  -  >  . 
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n  —  low  frequency  magnitude  of  noise  spectrum 
wQ  =  half-power  frequency  of  noise  spectrum 
If  the  frequency  output  of  the  Gaussian  random  generator  varied 
from  0  radians/sec  to  00  radians/sec,  the  value  of  w  would 
be  one.  Since  neglecting  the  high  frequency  filter  on  the 
generator  is  the  same  as  assuming  the  frequencies  go  to  00 , 
the  value  of  w^  is  taken  as  one.  (In  all  the  following  cal¬ 
culations  one  may  note  that  they  are  not  highly  sensitive  to 

changes  in  w  values  as  the  power  of  w  is  the  same  in  both 

o  o 

numerator  and  denominator.)  To  obtain  the  value  of  n  in  the 
Laplace  transform  representation  of  the  Gaussian  random  genera¬ 
tor  output,  consider  the  new  block  diagram  of  the  mathematical 
model  used  on  the  computer  (see  Figure  11) . 


I  reference 


Figure  11 
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On  the  computer  the  control  loop  was  broken  along  the 
dashed  line  in  Figure  11,  The  Gaussian  random  generator 
was  then  adjusted  until  the  rms  of  h  above  the  dashed  line 
was  one,  (The  Gaussian  random  generator  was  then  left  at 
this  setting  for  all  random  studies,)  With  the  open  loon  h 
equal  to  one,  "n"  is  calculated  as  follows: 

For  the  open  loop: 


nw 


h(s) 


o 


-)  ( 


k9 

•)  (— > 


K~nw 
2  o 


(18) 


S  +  w  t  S  +  1  S 
o 


T S  +  ( tw  +  1 ) S  +  w 
o  o 


ms 


By  equations  (10)  and  (11) 


(K0nw  ) ^  j  dw 
£  o 


(19) 


2uj  -«  {^(jw)  +(  w  +  l)jw  +  w  1  f t  ( - jw)  “+  ( tw  +.l)(-jw)+w  } 

J  o  o  o  o 


Using  the  integral  tables  (5-1)  of  reference  (9) 


h 


(K.nw  ) 
2  o 


ms 


2w  (tw  +  1) 
o  o 


(20) 


Since  the  results  are  going  to  be  compared  with  those  obtained 
on  the  computer  model,  the  constants K2  and  t  are  made  equal 
to  the  computer  values  in  the  analytical  calculations. 


Therefore , 


K2  =  0,2 


t  =  10 


w  -  1 

o 


Thus  by  equation  (20) 


h  =  1 

rms 


h 


ms 


(0c2n) 
2(10  +  1) 


1 


Therefore , 


n  -  23,4 


. 
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The  Laplace  transform  of  the  Gaussian  random  generator  output 

23 «  4 


is  thus 


(In  the  calculations  that  follow,  all  ex- 


S  +  1 

pressions  will  be  derived  before  values  are  substituted. 


Therefore,  the  Gaussian  random  generator  output  will  be 

nw 

represented  by  QR(S)  =  - —  until  values  are  again  substituted.) 


S  +  w 


Now  with  Q  (S)  as  the  Gauosian  random  generator  output, 
the  transfer  function  is  obtained  for  each  variable  in  the 
closed  control  loop  of  Figure  11.  All  average  values  m  the 
mathematical  model  are  taken  as  zero.  Therefore,  the  reference 


value  in  Figure  11  is  zero. 
Therefore , 

From  the  block  diagram  of  Figure  11 


E (S )  -  hf(S) 


K. 


hf(S)  = 


S (TpS  +  1) 


Op  (S )  s 

{  _JL -  -  Q(S)  } 

ts  +  1 


Q(S)  =  Ho  K  E(S) 


Also  , 

Now  substituting  equation  (21)  and  (23)  in  (22) 


(22) 


(23) 


E(S)  = 


Qr(S)  k2 


(tS  +  1)  (T_S  +  1) 

r 


Rearranging  equation  (24) 


K  KAl  E  (S) 

Z  3 _ 

(TpS  +  1)  S 


(24) 


E(S)  = 


k2s  Qr(S) 


(tS  +  1) (TpSZ  +  S  +  KxK2Ho) 


(25) 


nw 


Letting  0^(3) 

K 


S  +  w 


and  multiplying 


V  wo  >fi£  5  r  i:  J  v:  soft.-.:ddT 


w*i 


w  +  8 

.  il  /  •[ 


.oy.es  si  IX  saL/pi1?  ru  tuXsv 
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where 


E  (S)  = 


B 


o 


B. 


B. 


B. 


B 


By  equation  (23) 


K0nw  S 
2  o 


BqS4  +  B-jS3  +  B2S2  +  B3S  +  B4 


T  t 
F 


(T  +  T_)  +  W  T  T 

F  o  F 


1  +  t KK0H  +  w  (T„  +  i) 
2  c  o  F 

KK0H  +  (i  +  tK  K0H  )w 
2  a  2  o  o 


w  KK~H . 
O  2  6 


(26) 


Q  (S)  = 


KK0H  nw  S 
2  o  o 


B  S4  +  B. S 3  +  B0S2  +  B0S  +  B , 
0  12  3  4 


(27) 


Also,  from  the  block  diagram 

h  (S) 


hf(S)  = 


TpS  +  1 


Therefore,  by  equation  (21)  and  (26) 
h(S)  =  - 


K2nwQ(TFS  +  1)S 


B  S4  +  B,  S 3  +  B0S2  +  B-jS  +  3. 
o  1  2  3  4 


'  2  8 ) 


(29) 


Let  E  ,  Q  ,  and  h  equal  the  mean-square  of  the 
error,  manipulated  variable,  and  controlled  variable  res¬ 
pectively.  Then  apply  equation  (10)  to  equation  (26),  (27), 
and  (29)  to  obtain  the  above  three  variables  as  functions  of 
H^,  the  statistical  describing  function. 
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Solving  equation  (26)  for  E  and  thus  E 

ms  rms 


ms 


w 

/  {  (K2nw  )  2  ( jw)  (- jw)  dw)/{|jB  (jw)4  +  Bn ( jw) 

IT  -00  ^  u  o  1 


2 

+  B2(jw)  +  B3(jw)  +  B4  J&o(-jw)4  +  B1(-jw) 


+  B2(~jw)2  +  B3  ('-jw)  +  B4  } 


(29) 


From  table  (5-1)  of  reference  (9) 


E 


(K~nw  )  B, 
2  o  1 


^  2(B1B2B3  -  BoB32  -  B12V 


(30) 


Now  substituting  in  the  numerical  values  of: 

w  =  1,  n  =  23 „ 4 ,  K  =  4,  K0  =  0.2,  x  =  10,  T_  =  0.5 


E 


mS  1.019  +  8.3H  +  4.38H  2 

a  a 


(31) 


There fore , 


rms 


1.019  +  8.3H  +  4.38H  ‘ 

o  o 


(32) 


Keeping  all  other  values  the  same  but  changing  K  to  1.9 


ms  1.019  +  3.94H  +  0.932H  2 

a  a 


(33) 


Therefore , 


rms 


1.019  +  3  a  94H  +  0.932H  ‘ 

o  a 


(34) 


>  «  ‘  i  -'Ik  +  M).a  +  <>  t)-3  ■* 


<ftr  IU  -  (*0*  -  £H,  i«s- 
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Now  solving  equation  (27)  for  Q  and  Q 

^  ms  rms 


00 


Q 


ms 


2  tt  - 


/ 


{(KK2H  cnwQ) 2 ( jw)  (- jw) dw  }  /  {  |Bq ( jw) 4  +  B1 ( jw) 3 


+  B2(jw)2  +  B^  ( jw)  +  B4i  [Bo  (- jw)  4  +  B1  (- jw)  3 


.  v  2 


+  B2(-jw)  +  B3(-jw)  +  B4J  } 


(35) 


From  the  integral  table  (5-1)  of  reference  (•) 


Q 


ms 


(KK0H  nw  )  Bt 
_ 2 _ O _ 1 _ 

2<B1B2B3  "  BoB3S  -  B12V 


(36) 


Again  substituting  the  above  numerical  values  with  K  =  4 

2 


H 


•ms 


0.0635  +  0  o  519H  +  0.274H 


(37) 


Therefore , 


H 


a 


rms 


\J~0~. 


0635  +  0.519H  +  0.274H 


a 


(38) 


For  the  above  values  with  K  =  1.9 


H 


Q 


a 


ms 


0.2835  +  1.095H  +  0.2595H  ‘ 

a 


Therefore , 


(39) 


H 


Q 


a 


rms 


\f~ o7 


2835  +  1 o  095H  +  0.2595H 


(40) 
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Now  solving  equation  (29)  for  h  and  h 

ms  rms 


ms 


2* 


—  oc 


i  (K2nwQ)  2  (Tp ( jw)  4-  (jw)2)  dw>/{  LBq  ( jw)  4  + 


+  B2(jw)2  +  B 3  ( j  w  j  +  B4  [_B 0  ( —  j w )  4  +  3^  ( -  jw)  3 


+  B2 (-jw)  +  B^ ( - jw)  +  B4  ]  ) 


(41) 


From  the  integral  table  (5-1)  of  reference  (9 


h 


ms 


(K2nwQP  (TFB3  +  B±) 
2(B1B2B3  -  BqB32  -  Bi234) 


(42) 


Again  substituting  the  above  numerical  values  with  K  -  4 


15  »  75  +  2  „  2H 


h 


a 


ms 


15  c  8  +  12  8  „  9H  +  6  8.  OH  " 

a  o 


(43) 


Therefore , 


15,75  +  2  c  2H 


rms 


15.8  +  12  8, 9H  +6  8,011 

o  o 


(44) 


The  above  values  with  K  =  1.9 


15.75  +  1.046H 


h 


ms 


15. 8  +  61. 1H  +14, 49H  ‘ 

o  o 


(45) 


Therefore , 


15  e  75  +  1.046H 


h 


rms 


15.8  +  61. 1H  +  14.49H  ‘ 

a  o 


( 46 ) 
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Thus,  Erms/  Qms  an<^  ^rms  are  a-^  <?:*-ven  as  functions 
of  H  ,  the  statistical  describing  function.  The  statistical 
describing  function  must  now  be  determined  in  terms  of  the 
input  rms  for  the  various  nonlinearities.  It  can  be  seen  from 

Figure  11  that  this  input  rms  value  to  the  nonlinearity  is  E 

u  rms 

For  simplicity,  assume  E  to  be  a  biased  estimate  of  a.  Now 
r  J  rms 

since  the  average  value  of  e  is  zero,  the  rms  and  standard 

deviation  are  equal.  Therefore,  it  remains  to  calculate  f ( v) 

and  to  apply  equation  (17)  for  each  nonlinearity  to  obtain  the 

statistical  describing  functions. 

The  various  nonlinearities  for  which  H  ,  the  statis- 

cr 

i  i  3 

tical  describing  function,  will  be  obtained  are  e|e|,  e  , 

|e|e^,  and  | e |  .  The  different  H^'s  will  then  be  studied 

e 

for  their  effect  on  the  control  loop  of  Figure  11. 

Calculation  of  H  for  e|e| 

a 

Let  (e|e|)^  =  v^  to  comply  with  the  above  symbols  in 

the  theoretical  discussion. 

Applying  equation  (16) 

2  2 
4  -v  /2cr 
oo  v  e 

fTvT  =  J  - — —  dv 

— oo  a  /2  TT 

■5  4 

=  3a 

Applying  equation  (17) 
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Therefore , 


=  E  v'  3 
0  rms 


Calculation  of  H  for  e 

a 

Again  let  (e  )  =  v 


Applying  equation  (16) 


f  (v) 


00 


—  00 


c  2  2 
6  -v  /2a 
v  e  ' 


a  /2  Tt 


dv 


15a 


Applying  equation  (17) 


H 


15  6 
—  0 


=  2  o  74  o' 


Therefore , 


H  =  2.74  (E  ) 

o  rms 


Calculation  of  H  for  |e|e 


Again  let  (|e|e^)^ 


v 


8 


Applying  equation  (16) 


f  (v)  - 


2  2 
8  -v  /2c 

v  e 


o  /2n 


dv 


-  52.5a 


8 


Applying  equation  (17) 


H  =  7.2  o' 

a 


Therefore , 


H0  =  72  (Erms> 


(  47  ) 


(49> 
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Calculation  of  Ha  for 


e 


3/2 


Again  let  ( 


3/2 


v 


Applying  equation  (16) 


f  (v)  = 


v  e 


2  2 
-v  /2a 


dv 


/2tt 


a 

/2tt 

Applying  equation  (17) 

0.639 

/a 

0.639 

— -  (50) 

/~E 

rms 

The  solutions  of  equations  (47)  to  (50)  with  equations 

(32)  and  (34)  were  carried  out  graphically.  Also,  the 

equations ( 38 )  and  (40)  were  plotted  so  that  the  effects  of 

the  different  solutions  could  be  seen  on  the  manipulated 

variable,  Q  .  Since  h  and  E  curves  follow  one 
inn  s  irm  s  mn  s 

another  closely,  only  the  Erms  curve  was  plotted.  (The  impor¬ 
tant  points  on  the  h  curve  are  marked  by  a  cross  (+).) 

This  extension,  of  including  the  manipulated  variable  curve 
along  with  values  of  hms  and  the  two  curves  of  ErmQ  vs.  Hq 


H 


Therefore , 


H 


(  f )  noxifiL  p©  -  r  y  ■  ~A 


(VI)  ncj :■  _  i.'  '  qA 

(  £  '• 


-  I  oi  (V  J  Jrti  •  -.}  to  •  iox/i/ioe  sr! r 

<01  '«•  >  Aao  b-.riAfio  j  a  w  ( •••  Z )  for»3  (£* ) 

'  i1  i  ?  •  ■  ’ "  >  )  r.o.: 

to  n  3d  blxroo  s/:ojLit/Ioe  Lb  arid' 

.  Q  *  eXdxsi^fiv 

■ 

■  ■■*  ■  t  :  n  in  * 

•  >i  •!. i  i/cr  •  rifiilT  n.ubwion  r  '  O  %  flOX  ?  •*  i'ix©  sxriT 


'  :  s-  or'  1 


34 


needed  for  solution,  is  important 0  The  reason  is  that  by 
plotting  the  curves  of  different  nonlinearities,  one  can  pre¬ 
dict  the  effect  the  nonlinearity  will  have  on  both  Q  and 

J  yrms 

^rms°  gives  the  designer  the  power  to  choose  those  non- 

linearities  which  will  be  suitable  from  an  rms  standpoint. 

Further  elimination  can  then  be  carried  out  using  other  criteria, 
such  as,  lowest  IT AM. 


Each  nonlinearity  for  which  a  statistical  describing 

function  was  developed  above  will  now  be  discussed  in  terms 

of  i.;s  effect  on  h^  .  the  controlled  variable,  and  Q  ,  the 

rms  rms 

manipulated  variable.  First,  note  Figure  15  on  which  equation 

(32)  and  (44)  are  plotted  Figure  .15  illustrates  how  closely 

the  h  and  E  curves  follow  one  another.  Also,  on  Figures 
rms  rms 

12A,  12B,  13  and  14,  the  points  (b)  and  (d)  are  linear  values 
of  E and  Q ,  respectively  and  the  points  (a)  and  (c)  are 
the  values  of  E  and  Q  respectively  produced  by  the  non- 

linearity  under  discussion.  Also,  all  graphs  are  for  a  gain  of 
four  unless  otherwise  stated. 


The  nonlinearity 


which  is  described  by 


e 

equation  (50),  has  a  solution  at  point  (aj)  on  Figure  14.  Point 
(c)  on  Figure  14  is  the  corresponding  Qrms  value  that  would 
result  if  the  nonlinearity  were  used  in  the  control  loop. 
Thus,  as  compared  to  the  linear  system,  this  nonlinearity, 

le3/2| 


rms 


ctnd  h 
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,  brings  about  a  small  decrease  in  E 
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FIGURE  13 
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and  a  small  increase  in  Qrms«  This  is  unacceptable  from  the 
standpoint  of  this  study  and  thus  this  nonlinearity  was  not 
given  further  consideration. 

The  nonlinearity  elei  is  shown  in  Figures  12A  and  12B 

with  solutions  at  points  (a)  and  corresponding  Q  values  at 

r  rms 

points  (c) .  As  seen  on  both  Figures,  this  is  an  acceptable 

nonlinearity  because  a  decrease  in  Qms  is  brought  about  by  an 

increase  in  and  hrmg .  Figure  12A  is  for  the  control  loop 

operating  with  a  gain  of  4.  Figure  12B  is  for  a  gain  of  1,9. 

The  reason  for  the  two  different  gams  will  be  seen  when  the 

ITAM  values  for  the  different  nonlmearities  are  discussed. 

3 

The  nonlinearity  e  ,  which  is  described  by  equation  (48) , 

is  shown  on  Figure  13  with  a  solution  at  (a)  and  corresponding 

Qrmo  value  at  (O.  Again  this  nonlinearity  is  acceptable 

from  rms  viewpoint  for  E  and  h  are  increased  for  a  de- 

crease  in  Qrms •  Note  that  this  change  is  greater  than  for 

the  nonlinearity  e|ej  on  Figure  12A.  However,  the  author 

3 

will  show  later  that  e  and  higher  powers  are  not  acceptable, 
in  his  judgement,  because  of  the  resulting  ITAM  values  pro- 

i v'  i  )  i  l  .  T'  .  :  .] 

duced  by  the  high  power  functions.  Therefore,  the  Figure 

'il  8  3  0  0 

with  the  solution  for  jeje  is  omitted. 

As  mentioned  above,  the  h  and  E  curves  shown  on 

rms  rms 

Figure  15  follow  one  another  closely.  The  difference  between 
the  two,  however,  is  the  amount  the  rms  is  lowered  by  re¬ 
moving  the  high  frequency  load  disturbances. 
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The  above-mentioned  figures  will  be  discussed  again 
in  4-he  experimental  section. 

Summary 

^his  section  has  done  basically  three  things.  First, 
a  mathematical  model  for  the  system  was  produced.  Second, 
the  system  was  checked  for  stability  under  linear  conditions 
and  two  cases  of  nonlinear  conditions.  Finally,  a  method  of 
calculating  rms  values  for  linear  and  nonlinear  control  loops 
was  explained.  The  method  was  then  used  to  calculate  and 

;  If' 

predict  rms  value  changes?  due  to  giv#h  nonlinearities.  How¬ 
ever,  no  ITAM  values  were"-  calculated .  These  are  taken  with 
the  analogue  computer  in.  the '  next12  s'edtrbn . 
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Experimental  Report 

».  t  ■  % 

The  following  experimental  report  is  divided  into  two 
sections.  The  firgt..  section  is  concerned  with  the  mathematical 
model  programmed  for  the  analogue  computer.  With  the  computer 
model,  the  analytical  rms  calculations  are  checked  and  I TAM 
values  are  gathered  for  the, different  types  of  control.  The 
second  section  deals  with-  ^ork  done  to  verify  the  mathematical 
model  used  in  this  study.  Actual  field  equipment  was  set  up 
and  operated  for  this  part  of  the  stpdy.  The  two  reports  that 
follow  are  written  as  individual  laboratory . studies . 

A.  Computer  Studies 

Purpose  -  The  computer  studies  were  ..carried  out  toc 
check  the  mathematical  calculations  in  the  theory  ’  section- ^ 
and  to  calculate  ITAM  values  ^for  different  gain  and  step  , 
sizes.  From  the  data  gathered,  it  would  then  be -possible 
to  conclude  which  nonlinear  function,  if  any,  was  the  best 

■J&' 

for  the.  system  of  Figure  (4)  . 

Apparatus  -  The  equipment  used  was  a  75  amplifier  Pace 
analogue  computer  and  an  8  channel  Pace  recorder  and  a  Gaussian 
random  signal  generator  (type  RG77)  manufactured  by  Servomex 

« l  I 

Controls  Limited.  (The  frequency  band  of  tha  generator  is 

yv 

shown  on  Figure  16A. ) 

Procedure  -  The,  procedure  for  the  work  done  on  the 

,  ...  <i  ■■  1 1  ■'  4.  '  S- 

computer  was  divided  intg.  .three  operations:  one,  to  set  the 
program  up  generally;  two./  ,  to  set  up  and  opejpate  the  program 
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Bode  amplitude  plots  showing  frequency  band  widths 


RANDOM  GENERATOR  FREQUENCY  BAND 

FIGURE  16  A 


FREQUENCY  BAND  INPUT  TO  COMPUTER  MODEL 

FIGURE  16  B 


FREQUENCY  BAND  INPUT  TO  FIELD  EQUIPMENT 
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Computer  itam  values  vs  gain 

FIGURE  17 


/ 


2 


-i_ - 1 - -J - - — l- 

3  4  5  6 


GAIN  K 


\  \f\  ia  *  css 


I  T  AM  VOL  T5  x  0.5 


STEP  SIZE 


-  -15  - 


ItAM  VALUES  VS 

FIGURE  18 


STEP  SIZE  VOLTS 


/  TAM  volts 


46  - 


9,Or 


8.0 


7.0 


6.0 


i 


$ 

k 


5.0 


4.0 


50 


2.0 


1.0 


Field  equipment  itam  values 

vs  GAIN 


FIGURE  19 


Linear  Control 


O  Nonlinear  Control  elel 


VALUES  OBTAINED  ON  FIELD  EQUIPMENT 
ALL  STEP  SIZES  3  VOLTS 


i  i . . —i - 1 - 1 

2  3  4  5  6 


GAIN  K 


. 


w 


-  4  7  - 

LINEAR  SYSTEM  WITH  RANDOM  GENERATOR  S 
MEAN  SQUARE  CIRCUITS . 

FIGURE  20 
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for  ITAM  studios;  and  three,  to  set  up  and  operate  the  program 
for  rms  studies: 

1)  The  analogue  computer  programs  for  this  study  were 
scaled  to  be  ten  times  as  fast  as  the  actual  process.  This 
same  time  scaling  had  then  to  be  applied  to  the  Gaussian  random 
generator.  Thus,  the  input  which  was  fed  into  the  computer 
model  is  that  shown  in  Figure"  1GB  while  that  fed  to  the  field 

e<3uipment  is  shown  in  Figure  16C.  No  magnitude  scaling  was 
done. 

2)  To  determine  ITAM  values,  the  computer  programs 
were  as  shown  on  the  flow  charts  of  Figure..  21  for  the  linear 
system  and  of  Figure  23  for  the  nonlinearity  e  .  ITAM  values 
were  obtained  for  the  nonlinearity  e|e|  using, the  program 
shown  on  Figure  22  with  switches  (1)  and  (2)  open  and  (3)  and 
(4)  closed. 

With  the  computer  components  as  arranged  in  the 
above  flow  charts,  the  following  procedure  was  used  to  gather 
ITAM  data  for  Figures  17  and  18.  First,  the  step  size  was 
fixed  at  3  volts  for  the  linear  program.  (The  reason  why  3 
volts  was  used  will  be  explained  later.)  The'  gain  was  then 
changed  to  obtain  the  information  for  the  "lineiar  control" 
curve  shown  on  Figure  17.  This  was1  then  repeated  for  the 
two  nonlinearities  e|e|  and  e  .  These  results  are  also  plotted 
on  Figure  17.  The  lowest  and  thus  the  optimum  ITAM  value 
for  the  linear  control  loop  was  then  noted  on  Figure  17.  (It 
occurred  at  a  gain  of  four.)  The  smallest  gains  which 
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produced  the  same  ITAM  values  as  the  linear  optimum  for  the 
two  nonlinearities  e|e|  and  e  were  then  chosen  from  Figure  17. 
fhe  gain  value  of  1.9,  chosen  for  the  nonlinearity  e|e|,  was 
then  set  on  the  program  shown  on  Figure  22.  Similarly,  the 
gain  of  0o65,  chosen  for  the  nonlinearity  e  ,  was  set  on  the 
program  shown  on  Figure  23.  With  the  gains  fixed,  the  step 
sizes  for  both  the  programs  shown  on  Figure  22  and  23  were 
changed  to  gather  the  data  for  Figure  18.  This  completed  the 
ITAM  studies. 

3.  In  view  of  the  results  shown  on  Figures  17  and  18, 
the  decision  was  made  to  drop  the  nonlinearity  e  and  to  use 
the  nonlinearity  e|e|  with  a  gain  of  1.9.  Thus,  rms  values 
were  only  calculated  for  the  linear  system  and  the  nonlinear 
system  with  e|e|.  This  was  done,  on  the  computer,  using  the 
components  as  arranged  in  the  flow  charts  of  Figure  20  and  22 
respectively.  For  rms  calculations  switches  (3)  and  (4) 
were  open  and  (1)  and  (2)  were  closed  on  Figure  22. 

For  the  rms  calculations,  the  Gaussian  random 
generator  had  to  first  be  set  so  as  to  give  a  value  of  n,  as 
mentioned  in  the  theory  section,  of  23.4.  This  was  done  by 
using  the  open  loop  circuit  that  is  contained  within  the 
dashed  rectangle  on  Figure  20.  Using  only  the  above  open 
loop  circuit,  the  Gaussian  random  generator  amplitude  was 
adjusted  until  the  rms  value  of  h  was  1.  By  the  calculations 
in  the  theory  section,  n  was  then  23.4.  With  the  Gaussian 
random  generator  left  at  this  setting,  the  control  loop  of 


>  3  p.b  dsufev  MAT  I  emsa  arlj  bsouboi 

. 


1  '  •'  ii  er  Sipcti  ■  >  h>,'i  .  V  . 


[  3  /pi  !  i '  '  *1  s  L*/yp  oa  >1 pn&rto 


^nt  -^ip-'  n  n  axi  sd..  boisiuolBO  vl  '  <  isw 


nfii  e  sr’i  anoi^slijolj  o  r.Jtn  £  w  xcl  ■ 

is  95  b9cj  ed  i2*j  ■]  o  -  &n  loi  > tonsp 


"»n  ?no3  ei  J-jBili  iic/o^io  qooi  n^jqo  ©rid  pni&u 

•  0£  Oi  '  i  .  9^  6--«  ' 

■ 


■ 


. 


5  3 


Figure  20  was  then  closed  and  the  rms  values  throughout  the 
loop  were  measured.  This  same  procedure  was  then  repeated 
for  the  nonlinear  system  shown  on  Figure  22.  This  ended  the 
rms  investigations. 

Results  -  The  ITAM  values,  as  mentioned  above,  were 
plotted  on  Figures  17  and  18.  The  fluctuations  of  the  mani¬ 
pulated  variable  and  controlled  variable  from  which  the  rms 
values  were  calculated  are  shown  on  Figures  28  and  29,  respec¬ 
tively.  Figure  28  compares  the  movement  of  the  manipulated 
variables  for  the  linear  system  and  the  nonlinear  system  e  ej 
under  the  same  load  disturbance.  Figure  29  compares  the  con¬ 
trolled  variables  f  r  • he  same  systems  under  the  same  load 
disturbance.  The  differences  between  the  two  types  of 
control  were  such  that  the  nonlinear  system  had  a  manipulated 
variable  rms  39%  lower  than  the  linear  system  and  a  cor  trolled 
variable  rms  67%  higher. 

Conclusions  and  Discussion  -  The  ITAM  results  showed 

that  the  nonlinear  functions  elej  and  e  were  first  worse 

than  the  linear  controller  (see  Figure  18)  then  better  and  then 

worse  again  as  the  step  change  sizes  were  increased.  The 

3 

results  also  showed  that,  e  varied  more  from  the  linear 
controller  than  e|ej  by  a  considerable  amount.  The  non- 
linearity  e'  had  a  shorter  range  of  step  sizes  than  eje  for 
which  its  ITAM  values  were  lower  and  then  better  than  the 
linear  controller  values.  Thus,  because  of  the  desire  to 
have  low  values  of  ITAM  for  a  wide  range  of  step  sizes,  e I e ; 
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was  chosen  and  e  and  higher  powers  discarded.  The  decision 
to  discard  the  higher  powers  was  also  influenced  by  the  fact 
that  the  I TAM  values  increased  very  rapidly  above  6  volts 
and  just  below  3  volts  as  seen  on  Figure  18. 

The  rms  values  of  the  linear  and  nonlinear  (eje|) 
showed  that  the  analytical  predictions  were  quite  good.  A 
38%  decrease  in  the  manipulated  variable  was  predicted 
(difference  between  (d)  on  Figure  12A  and  (c)  on  Figure  12B) 
for  an  increase  of  64.5%  in  the  controlled  variable  (differ¬ 
ence  between  cross  (+)  on  Figure  12A  and  cross  on  Figure  12B) . 
Instead  a  39%  decrease  was  obtained  in  the  manipulated 
variable  for  an  increase  of  67%  in  the  controlled  variable. 

The  actual  voltage  values  of  the  manipulated  variable  and 
controlled  variable  agreed  within  5%  of  the  analytical  values. 
Both  values  were  low  on  the  computer.  This  could  be  due 
partly  to  the  fact  that  the  high  frequency  filter  in  the 
Gaussian  random  generator  was  neglected  in  the  analytical 
calculations . 

The  computer  results  show  that  e|e|  is  a  successful 
nonlinear  controller.  If  the  nonlinear  controller,  e|e|,  is 
set  properly,  it  can  give  a  lower  ITAM  value  than  the  linear 
controller  and  still  lower  the  rms  of  the  manipulated  variable 
considerably ,  Also,  should  the  step  size  be  outside  the  region 
where  the  nonlinear  controller  ITAM  value  is  higher  than  that 
for  the  linear  controller,  the  difference  will  not  be  drastic. 
All  that  remains  now  is  to  confirm  the  above  results  on  the 
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B„  Field  Equipment  Studies 

Purpose  -  The  purpose  of  conducting  the  tests  on  the 
field  equipment  was  to  check  the  accuracy  of  the  mathematical 
model  and  to  see  if  the  operations  as  described  on  the  com¬ 
puter  could  be  carried  out  on  the  field  equipment. 

Apparatus  -  The  equipment  used  was  as  follows:  an 
Applied  Dynamics  32  amplifier  analogue  computer,  two  Foxboro 
needle  valves  with  pneumatic  motors,  one  Foxboro  differential 
pressure  cell  for  liquid  level  measurements ,  a  Gaussian  random 
signal  generator  (Type  RG77)  made  by  Seromex  Controls  Limited, 
and  three  Honeywell  current  to  pressure  and  pressure  to 
current  transducers  for  communication  between  the  process  and 
the  computer. 

Procedure  -  The  procedure  again  consisted  of  three 

steps . 

1)  The  first  step  was  to  assemble  the  apparatus 
such  that  any  one  of  the  arrangements  shown  on  Figure  24, 

25,  26  and  27  could  be  coupled  together  easily.  With  the 
field  equipment  there  was  no  time  or  magnitude  scaling.  Be¬ 
cause  of  this,  the  time  constant  on  the  low  frequency  filter 
was  reduced  by  a  factor  ten  from  the  time  scalec  computer 
value  (see  Figure  1.6C)  ,  This  was  done  to  make  the  inputs  to 
the  computer  model  and  the  field  equipment  approximately  tne 
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ON  FIELD  EQUIPMENT 
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2)  To  obtain  the  ITAM  values  the  computer  and  pneumatic 
equipment  were  set  up  as  shown  schematically  in  Figures  24  and 
26. 

For  the  linear  system  of  Figure  24,  the  potentio¬ 
meters  were  first  fixed  so  that  the  tank  was  half  full  of 
water.  The  flow  of  water  through  the  system  was  then  set  at 
50  percent  of  its  maximum  value „  (This  was  done  by  adjusting 
the  potentiometers  marked  qA  on  Figure  24.)  The  reset  potentio¬ 
meter  was  then  adjusted  until  zero  error  was  obtained.  A  3- 
volt  step  was  then  set.  Finally,  the  3-volt  step  was  applied 
and  the  potentiometer  turned  on  that  is  set  at  0.01.  The  ITAM 
values  were  then  read  as  shown  on  Figure  24.  This  same  pro¬ 
cedure  was  repeated  for  different  gains  to  collect  the  data 
for  Figure  19. 

For  the  nonlinear  system,  eje|,  of  Figure  26,  the 
procedure  for  obtaining  ITAM  data  was  the  same.  The  results 
collected  for  the  nonlinear  system  are  also  shown  on  Figure  19. 

3)  For  obtaining  rms  values,  the  equipment  was  set 
up  as  shown  on  Figures  25  and  27, 

The  linear  system  (see  Figure  25)  again  had  its 
tank  level  and  water  flow  adjusted  so  that  they  were  50 
percent  of  their  maximum  values.  Next  the  control  loop  was 
opened  at  the  input  to  the  feedback  filter  where  the  potentio- 
ijigter  is  marked  0.2  on  Figure  25.  With  this  open,  loop  circuit, 
the  Gaussian  random  generator  was  adjusted  until  an  hrms 
one  was  obtained.  The  gain  was  then  set  at  2,0.  (This  was 
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Figure  29. 
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the  gain  that  produced  the  lowest  I TAM  value  for  the  linear 
system  on  Figure  19 0)  The  control  loop  was  connected  again  and 
the  response  of  the  controlled  variable,  h,  and  the  manipulated 
variable,  q  ,  recorded. 

The  above  procedure  was  repeated  again  for  the  non¬ 
linear  system,  e|e| ,  with  the  components  as  shown  on  Figure  27. 

The  only  difference  was  that  a  gain  of  1.05  was  used.  This 
was  so  the  ITAM  values  matched,  as  on  the  computer  model,  at  a 
step  size  of  3  volts. 

Results  -  The  ITAM  values  are  shown  on  Figure  19. 

The  recordings  of  the  controlled  variable,  h,  and  the  mani¬ 
pulated  variable,  q  ,  are  shown  on  Figures  30,  31,  32  and  33. 
Figure  30  shows  a  comparison  of  qQ  as  it  was  measured  electrically 
before  the  transducer.  (Note  here,  however,  that  all  com¬ 
parisons  were  not  necessarily  under  the  same  load  disturbance 
as  the  linear  and  nonlinear  systems  could  not  be  run  simul¬ 
taneously  on  the  field  equipment  as  on  the  computer.) 

Figure  31  shows  a  comparison  of  h  for  the  linear  and  non¬ 
linear  systems  under  the  Gaussian  random  load  fluctuations. 

Figure  32  and  33  taken  together  compare  the  linear  and  non¬ 
linear  values  respectively  of  the  manipulated  variable,  q  ,  as 
measured  pneumatically.  The  rms  decrease  in  the  manipulated 
variable,  q  ,  as  measured  pneumatically  was  approximately  35 
percent.  The  increase  in  the  controlled  variable,  h,  was 
approximately  61  percent.  The  word  approximate  was  used 
because  the  results  change  slightly  if  different  sections  of 
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the  two  runs  with  the  linear  and  nonlinear  systems  are  compared. 

Conclusions  and  Discussion  -  The  field  equipment  studies 
showed  the  same  trends  and  approximate  results  as  the  theoret¬ 
ical  work. 

First,  note  how  the  ITAM  curves  from  the  computer  studies 
and  field  equipment  studies  show  essentially  the  same  relation¬ 
ship  between  the  linear  curve  and  the  nonlinear,  e|e|,  curve 
(see  Figures  17  and  19) .  In  both  cases  the  curve  for  the  non¬ 
linearity,  e|e|,  as  compared  to  the  linear  curve,  is  of  a 
higher  slope  and  reaches  its  lowest  value  at  a  smaller  gain. 
However,  one  can  also  see  that  the  actual  values  on  the  corres¬ 
ponding  curves  from  the  two  different  studies  do  not  match  in 
slope  or  gain.  The  nonlinear  valve  on  the  field  equipment  is 
the  most  likely  reason  why  the  above  discrepancies  occurred 
between  the  comparable  curves  on  Figures  17  and  19. 

Second,  11.5  percent  was  the  largest  percent  error  made 
in  predicting  rms  changes.  This  occurred  when  the  computer 
values  predicted  a  39  percent  change  in  the  manipulated  variable, 
q  ,  and  the  field  equipment  only  produced  a  35  percent  change. 

The  analytical  results  had  predicted  a  38  percent  change  and 
were  thus  out  by  a  9  percent  error.  This  was  the  greatest 
deviation,  also,  of  the  analytical  calculations  from  the  field 
equipment  results. 

Although  the  fact  that  the  suspected  nonlinearity  in  the 
valve  not  taken  into  account  in  the  theoretical  study  may  have 
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contributed  in  large  measure  to  the  discrepancy  between  theo¬ 
retical  and  experimental  results  and  consequently  detracted 
somewhat  from  the  effectiveness  of  the  mathematical  model  used. 
It  must  be  pointed  out  that  the  theoretical  study  using  the 
linear  model  did  serve  to  narrow  down  the  choice  of  controller 
for  the  system.  To  this  extent  then  the  model  can  be  con¬ 
sidered  to  be  useful. 
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General  Conclusions  and  Discussion 

In  this  final  section,  the  following  topics  are 
coveredc  First,  a  defense  is  given  of  the  three  criteria  that 
guided  the  outcome  of  this  study.  Second,  the  reasons  are 
presented  for  using  a  3-volt  step  size  to  compare  ITAM  values 
and  a  Gaussian  random  generator  setting  that  resulted  in  an 
open  loop  rms  value  of  h  equal  to  one  for  comparing  all  rms 
values.  Third,  problems  with  the  experimental  equipment  are 
discussed.  Fourth,  the  different  types  of  control  are  re¬ 
examined.  Finally,  the  author  states  why,  in  his  judgement, 
the  nonlinearity  e[e|  is  the  best  to  use  in  the  control  loop 
of  Figure  7. 

The  first  criteria  of  this  study,  stability, is 
necessary  for  any  controller.  Without  stability  there  is  no 
control  and,  therefore,  it  is  the  most  important  criteria  for 
any  controller.  The  second  criteria,  ITAM,  is  important 
because  industry  wants  start-up  disturbances  to  damp  out 
rapidly.  This  is  because  start-up  disturbances  are  generally 
the  worst  a  system  has  to  undergo  and  industry  does  not  want 
them  to  last  long  because  of  the  problems  they  create  with 
the  rest  of  the  down-stream  processes.  Since  the  manipulated 
variable,  qQ,  was  the  out-flowing  stream  in  this  study,  ITAM 
describes  the  above  desire  mathematically.  This  criteria  was 
placed  second  in  preference  because  of  the  importance  of  start¬ 
up  problems  in  Industry.  Industry  also  desires  the  controller 
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to  move  the  manipulated  variable  only  when  necessary  during 
actual  operation.  Again  this  is  so  that  the  downstream  pro¬ 
cesses  will  be  disturbed  as  little  as  possible.  This  desire 
is  represented  mathematically  be  demanding  a  low  rms,  (The 
average  flow  value  is  considered  as  zero  and  fluctuations  from 
zero  are  used  to  calculate  the  rms  value.) 

A  value  of  3  volts  as  the  comparison  voltage  for  all 
I TAM  values  was  chosen  because  the  author  felt  that  most  start¬ 
up  conditions  would  be  such  as  to  make  the  head  change  at 
start-up  three  times  that  of  the  open  loop  rms  value  of  the 
head  variation.  This  selection  was  quite  arbitrary,  however, 
and  others  may  prefer  some  other  value.  The  reason  for  using 
an  open  loop  rms  value  of  h  equal  to  one  was  that  this  pro¬ 
vided  a  definite  standard.  If  a  closed  loop  value  were  used 
the  type  of  controller  within  the  control  loop  would  have 
played  a  part.  The  reason  a  given  output  from  the  Gaussian 
random  generator  was  not  used  as  a  standard  was  because  of 
the  nonlinear  valves.  The  author  wanted  to  make  certain  the 
field  equipment  was  undergoing  approximately  the  same  sort  of 
load  disturbance  as  the  computer  model.  (Using  the  open  loop 
value  also  aided  the  analytical  calculation  in  that  "n"  for 
the  Gaussian  random  generator  was  then  easily  calculated.) 

The  author  had  to  tolerate  two  unknowns  in  this 
study  because  of  the  equipment  used.  First,  the  specifications 
supplied  by  the  manufacturers  of  the  Gaussian  random  generator 
indicated  that  the  frequency  band  below  0.25  radions/sec  was 
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uncertain c  The  manufacturers  said  they  believed  the  band  was 
flat  to  nearly  0C0  radions/sec«  but  they  could  not  prove  it. 

This  is  why  the  diagrams  of  Figures  16A,  16B  and  16C  are 
drawn  with  dashed  lines  below  0*25  radions/sec.  Also,  note 
on  Figures  16B  and  16C  that  the  band  widths  from  the  Gaussian 
random  generator  fed  to  the  computer  and  field  equipment  were 
not  equal.  Thus,  the  field  equipment  should  suffer  more  from 
high  frequencies.  This  was  born  out  by  the  data  of  Figure  31. 
The  second  unknown  the  author  had  to  tolerate  was  the  fact  that 
the  linear  and  nonlinear  systems  could  not  be  run  simultaneously 
on  the  field  equipment.  Thus,  one  could  not  be  certain  that 
the  load  disturbances  on  both  systems  were  alike.  However, 
because  fairly  long  runs  were  used  to  calculate  the  rms  values, 
they  are  probably  quite  precise. 

The  nonlinearity  e|e|  was  chosen  in  this  study,  over 

3 

higher  powers  such  as  e  ,  for  two  reasons*  One  was  because 
ejej  has  a  longer  range,  than  the  higher  powers,  in  which  its 
ITAM  value  was  better  than  the  linear  controller.  Two,  when 
the  ITAM  values  of  e|e|  were  greater  than  the  linear  con¬ 
troller,  the  difference  was  not  as  great  as  that  of  the 
higher  powers.  Thus,  the  gain  for  the  nonlinear  controller, 
e|e|,  was  not  as  difficult  to  set  to  be  in  the  valley  just 
above  the  3-volt,  step  size  on  the  curves  shown  on  Figure  18. 
Also,  as  can  be  seen  by  the  analytical  calculations,  the 
biggest  change  in  rms  is  going  from  the  linear  system  to 
the  nonlinear  system  with  e|e|.  When  one  goes  from  e|ej  to 
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e  ,  the  change  in  the  qQ  rms  is  smaller  than  above  and  for  the 
same  change  in  the  qQ  rms  value  a  greater  percent  increase  in 
rms  rec3uired<-  The  major  disadvantage  of  the  controller 
with  the  nonlinearity  e|ej  was  that  its  maximum  h  value  was 
35  percent  higher  than  the  linear  controller’s  maximum  h  value a 
(This  figure  was  obtained  by  considering  all  runs.)  Because 
of  this  one  might  argue  that  the  same  effect  as  that  given  by 
the  nonlinear  controller,  ejej,  could  be  obtained  by  lowering 

the  gain  on  the  linear  controller  so  that  its  h  was  35 

rms 

percent  higher.  Doing  this,  however,  does  not  lower  qQ  rms 
the  approximate  39  percent  obtained  with  the  nonlinear  con¬ 
troller.  Instead  only  an  approximate  30  percent  lowering  of 

the  rms  of  q  is  obtained  and  the  ITAM  value  would  be  increased 

o 

an  approximate  110  percent.  Thus,  the  start-up  response  would 
be  greatly  decreased  in  desirability  and  the  manipulated 
variable  would  not  be  lowered  as  much  as  that  by  the  non¬ 
linear  controller. 

The  author  feels  this  study  has  contributed  informa¬ 
tion  in  two  aspects.  First,  it  has  designed  a  nonlinear  con¬ 
troller  that  works  well  for  the  system  described  in  this 
study.  Also,  information  on  stability,  ITAM  values  and  rms 
values  have  been  given  for  the  controlled  system.  Second, 
the  study  has  shown  how  to  use  the  statistical  describing 
function  method  to  calculate  rms  values.  It  has  also  extended 
the  above  method  to  include  the  manipulated  variable,  qQ , 
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The  study  verified  the  statistical  describing  function  approach 
for  analytical  calculations  to  within  4e5  percent  error  on 
the  computer  model  and  to  within  9  percent  error  on  the  field 
equipment  0 

In  conclusion,  the  author  would  like  to  point  out  that 
eje|  was  chosen  as  the  nonlinearity  to  use  in  this  study 
because  it  fit  the  special  specifications  the  besto  However, 
should  the  demands  on  the  controller  change,  there  is  little 


likelihood  that  e  e 


would  be  the  answer  as  each  special 


problem  in  nonlinear  control  theory  demands  a  different  type 


of  function. 
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